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I. INTRODUCTION 



In physics of strong nuclear interactions, the observed hadrons at low energies can be 
described in terms of meson and baryon degrees of freedom, while at high energies the de- 
scription is formulated in terms of quarks and gluons. Quark-hadron duality is an intriguing 
phenomenon which, by connecting the two energy regimes, provides a dual description of 
hadronic variables. Its understanding, however, still represents a fundamental challenge. 

In Ref. jl| we present a novel new approach to study quark-hadron duality using truncated 
moments of structure functions, or the integrals of structure functions over restricted regions 
of Bjorken variable x. In particular, by studying the evolution of structure functions 
integrated over specific nucleon resonance regions, where is the virtuality of the exchanged 
photon, we determine the degree to which individual resonance regions are dominated by 
leading twist, and hence quantify their higher twist content. 

The phenomenon of quark-hadron duality was first observed in inclusive electron-nucleon 
scattering at SLAC by Bloom and Oilman in 1970 It refiects the similarity between 
structure functions averaged over the resonance region and the scaling, or leading twist 
function. The former is characterized by hadronic bound states, whereas the latter describes 
the high energy deep inelastic continuum, in other words, scattering from free quarks j^. 
The oscillation of nucleon resonances around the scaling curve shows that hadrons follow 
the QCD scaling behavior. The only rigorous interpretation of the so-called Bloom-Gilman 
duality within the theoretical framework has been done in terms of the moments of structure 
functions. They were analyzed within the QCD operator product (or twist) expansion. Here 
the leading term in expansion is given by matrix elements of leading twist quark-gluon 
bilocal operators, and is associated with free quark scattering. Furthermore, higher terms 
in the expansion correspond to nonperturbative (higher twist) quark-gluon interactions. In 
this language duality is explained as the suppression of higher twist contributions to the 
moments 

Experimental data from recent years ^ suggest that quark-hadron duality not only exists 
over the whole resonance region but also locally. It holds over restricted regions of hadronic 
final state mass W, or even for individual resonances. Nevertheless, any insight about 
workings of local duality has been confined to QCD-inspired models of the nucleon. As such 
our understanding of quark-hadron duality in nucleon structure functions within QCD is 



2 



incomplete. 

Many earlier analyses of duality within a QCD context have used full moments of structure 
functions and hence quantified the higher twist content over all region of x. Our approach, 
on the other hand, allows one to study for the first time the distribution of higher twist 
corrections over various regions in x (or W) in a well-defined systematic way. 



II. TRUNCATED MOMENTS AND EVOLUTION 

In analogy with full moments, e.g. of a parton distribution function (PDF) g(x, Q^), 

M„iQ^) = f'dxx^-' g(x,Q2) , (1) 
truncated moments are defined as: 







Mnixo,l,Q') = [ dxx^-' q{x,Q^) , (2) 

where the integration over the Bjorken x variable is restricted to Xq < x < 1. By truncating 
the integration region to some minimum value Xq, one avoids the problem of extrapolating 
parton distributions into unmeasured regions at small x. Several years ago this original 
idea was introduced to study structure function moments for which small-x data were not 
available Note that x is related to the invariant mass squared W'^ of the virtual photon- 
hadron system through W"^ = + (1 — x) /x, where M is the nucleon mass. Thus an 
infinite energy is required to reach the x ^ limit, and accordingly some extrapolation to 
X = is always needed in practice to evaluate the moment. Unlike full moments, truncated 
moments satisfy the evolution equations, which mix the lower moments with the higher ones. 

An alternative approach, the so-called diagonal formulation of truncated moments, was 
later developed to avoid the complication of mixing j?]. It was shown that truncated mo- 
ments obey the DGLAP evolution with a modified splitting function in the Mellin convolu- 
tion. Namely, 



dt 2-K 



iP:,0Mn)ixo,Q') , (3) 



where 



P:,{z,as{Q')) = z^P{z,as{Q')) (4) 
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plays the role of the splitting function for the truncated moments. Moreover, by introducing 
the doubly-truncated moment of g(a;, Q^), 

A4n(a;mm,a;max,Q^) = / dx x"^'^ q{x, Q"^) , (5) 

^min 

the evolution equations for truncated moments can be generalized to any subset in the 

X-region, X^in <X< Xmax- 

In our work we partially follow the latter approach and apply it to the study of the proton 
structure function F2 in the large-x region, populated by nucleon resonances. We use recent 
high-precision data from Jefferson Lab experiment E91-110 [8|, and quantify the size of the 
higher twists for the lowest three moments in various regions oi W < 2 GeV. By dividing 
the data into the three traditional resonance regions we extract the leading and higher twist 
content of each region. 



III. DATA ANALYSIS AND EXTRACTION OF HIGHER TWISTS 

In order to determine the degree to which nucleon structure function data in specific 
regions in x (or W) are dominated by leading twist we first evolve the structure functions by 
brute-force using a suitable numerical integration routine, apply the target mass corrections 
(TMCs), and finally calculate their corresponding truncated moments over the integration 
range Wth <W< PKnax, where PVth = M + is the inelastic threshold. Deviations of the 
evolved moments, computed to next-to-leading order (NLO) accuracy, from the experimental 
data at the new then reveal any higher twist contributions in the original data. In 
particular, we analyze recent data on from Jefferson Lab covering a range in from 
< 1 GeV^ to ^ 6 GeV^. 

The evolution of the measured truncated moments requires the structure function to be 
decomposed into its nonsinglet and singlet components, which are a priori unknown. We 
shall assume that in our region of interest (i.e. at moderate to large x) the structure function 
F2 is well approximated by its nonsinglet component, and hence will evolve the truncated 
moments as nonsinglets. To test the accuracy of the nonsinglet (NS) evolution versus the 
full or exact evolution we use the trial function, whose decomposition into its nonsinglet 
and singlet components is known. In our case this is the leading twist proton F2 structure 
function computed from the MRST2004 PDF fit {l^. First we evolve its nonsinglet and 
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Figure 1: Ratio of the truncated moments of F2 evolved from = 25 to 1 GeV^, using NS and 
full evolution, versus the truncation point Wmax for the n = 2 (solid), 4 (dotted) and 6 (dashed) 
moments. 

singlet components separately and then construct the full function, and secondly we evolve 
the trial function under the assumption that the total function can be treated as a nonsinglet. 
The comparison between the truncated moments calculated from the two evolutions gives 
the accuracy of the NS evolution of the various moments with respect to the exact results 
using both singlet and nonsinglet evolution. The ratios of these are plotted in Fig. [T] as a 
function of the truncation point W^max = + Q'^i^/^min — !)• In the traditional nucleon 
resonance region, 1.2 ^ W^max ~ 2 GeV, the differences between the full and NS evolution are 
of the order 2-4% and increase with increasing PVmax- The n = 2 moment is most sensitive 
to singlet evolution. As expected, the differences are smaller for the higher moments. Thus 
in the region relevant for our study one can therefore conclude that the evolution error, or 
the uncertainty introduced by evolving the truncated moments as nonsinglets, is less 
than 4%, which will be included in the errors in our final results. 

In the analysis we assume that the data beyond some large value of (in our 
case = 25 GeV^) are dominated by twist-2 contributions, which is consistent with most 
global analyses of PDFs. For example, a comparison of the DIS data fit shows the excellent 
agreement between the leading twist structure function, after correcting for TMCs, and the 
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Figure 2: dependence of the ratio of n = 2 truncated moments M-2 calculated from the data 
and from leading twist evolution from Qq = 25 GeV^ (including TMCs), for various intervals in 
W: the first (A) resonance region (diamonds), second (-Sii) resonance region (squares), the first 
and second combined, corresponding to W^^-^ = 2.5 GeV^ (dotted curve), third (-F15) resonance 
region (open circles), first three regions combined, W^max = 3.1 GeV^ (dashed curve), and the entire 
resonance region W^g^^ = ^ GeV^ (filled circles). Note that some of the points are offset slightly for 
clarity. 

data at this scale. Moreover, we observe that the inclusion of TMCs has to be properly 
accounted for before drav^^ing any conclusions about higher twists from data. To quantify 
the higher twist content of the specific resonance regions, and at different values of Q^, we 
consider several intervals in W: PV^h — < 1.9 GeV^, corresponding to the traditional 
A(1232) (or first) resonance region; 1.9 < < 2.5 GeV^ for the S'ii(1535) (or second) 
resonance region; and 2.5 < < 3.1 GeV^ for the Fi5(1680) (or third) resonance region. 
The n = 2 truncated moments corresponding to these regions are plotted in Fig. [2] for various 

values. It's worth noting that below = 1 GeV^ the applicability of a pQCD analysis 
becomes doubtful and the decomposition into leading and higher twists is no longer reliable. 

The entire resonance region data (filled circles in Fig. [2]) deviate from leading twist 
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Figure 3: dependence of the fractional higher twist (HT) contribution to the n = 2 truncated 
moment data, for various intervals in W (as in Fig. [2]). 

behavior at the level of ^ 15% for all values of with a significant dependence for 
^ 3 GeV^, which is made more explicit in Fig. [3l Here the higher twist contributions 
to J^2, defined as the difference between the total and leading twist moments, are shown 
as ratios of the moments evaluated from the data. The strong dependence of the higher 
twists is evident in the change of sign around = 2 GeV^, with the higher twists going from 
^ -10% at = 1 GeV^ to 10-15% for ~ 5 GeV^. The higher twists, as expected, 
decrease at larger once the leading twist component of the moments begins to dominate. 

Looking at the individual resonance regions we see that in the A region (diamonds) the 
higher twist contributions are smallest in magnitude at large Q^, while they are largest for the 
Sii region (squares). Combined, the higher twist contribution from the first two resonance 
regions (dotted curve) is ^ 15% in magnitude for all . Furthermore, the higher twist 
content of the F15 region (open circles) is similar to the Sn at low , however, it decreases 
more rapidly for > 3 GeV^. The higher twist content of the first three resonance regions 
combined (dashed curve) is ^ 15-20% in magnitude for < 6 GeV^. Finally, integrating 
up to W^max — 4 GeV^ (filled circles), the data on the n = 2 truncated moment are found to 
be leading twist dominated at the level of 85-90% over the entire range. 
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Figure 4: (a) dependence of the ratio of truncated moments M./^ calculated from the data and 
from leading twist evolution from Q\ = 25 GeV^ (including TMCs), for various intervals in W 
(labels as in Fig. [2]). (b) Fractional higher twist contribution to the n = 4 truncated moment data, 
for various intervals in W (as in Fig. [3]). 



The results in Figs. [2] and [3] contain the experimental uncertainty on the F2 data (sta- 
tistical and systematic), and the uncertainty from the NS evolution of the data. For the 
experimental error we take an overall uncertainty of 2% for all truncated moment data, 
except for the n = 4 and n = 6 moments at W^max = 1-9 and 4 GeV^. Here the experimental 
uncertainties are 4% and 3% for A^4, and 5% and 4.5% for M.^. The evolution error, on the 
other hand, is estimated by comparing the NS evolution with the full evolution using the 
MRST2004 fit, as in Fig. [H with the appropriate correction factor applied at each and 
W interval. 

For the higher moments, illustrated in Figs. [4] and [5], the overall magnitude of the higher 
twists is qualitatively similar to the n = 2 moments, although the values at which they 
start decreasing in importance are larger. In addition, the higher twist contributions at low 

are relatively larger for higher moments. For example, at = 1 GeV^ the magnitude 
of the higher twist component of the W"^ < 4 GeV^ region increases from ~ 10% for the 
n = 2 moment, to ~ 15-20% for n = 4, and ~ 20-30% for n = 6. This behavior can be 
understood from the relatively greater role played by the nucleon resonances in the large-x 
region, which is emphasized more by the (x-weighted) higher moments. 
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Figure 5: (a) Ratio of truncated moments M.^ calculated from the data and from leading twist 
evolution from Q\ = 25 GeV^ (including TMCs), for various intervals in W (labels as in Fig. [2|). 
(b) Fractional higher twist contribution to the n = 6 truncated moment data, for various intervals 
in W (as in Fig. [3]). 

IV. SUMMARY AND OUTLOOK 

Truncated moments of structure functions provide a firm foundation for the quantitative 
study of local quark-hadron duality within a perturbative QCD context. In this analysis we 
have precisely determined the degree to which duality holds as a function of in various 
resonance regions (i.e. the degree to which individual resonance regions are dominated by 
leading twist), and have calculated for the first time their higher twist content. In case of the 
n = 2 truncated moment of the proton structure function F2 we find that deviations from 
leading twist behavior in the resonance region iy^ < 2 GeV) are at the level of ^ 15% for 

> 1 GeV^. Morevover, there is a significant dependence in the ratio of moments of data 
to leading twist for ^ 3 GeV^, with the higher twists changing sign around = 2 GeV^. 
Our results indicate that in the second (5*11) and third (-F15) resonance regions the higher 
twist component is larger in magnitude than in the A resonance region. Similar behavior 
to the n = 2 truncated moments is found for the n = 4 and n = 6 truncated moments, 
however, here due to the large-x enhancement the resonances play relatively greater role, 
which then leads to larger higher twists at the same Q^. In addition, we have also quantified 
the uncertainty associated with evolving the structure function data as a nonsinglet, which 
was found to be ^ 4%. 
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The illustrated analysis is an encouraging new approach to understanding local Bloom- 
Gilman duality within a well-defined theoretical framework. It opens the way to further 
study of local duality in other structure functions, such as the longitudinal structure function 
Fl or spin-dependent structure functions. 
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